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ABSTRACT 

Mass-loss rate, M, is one of the key parameters affecting evolution and observational 
manifestations of massive stars, and their impact on the ambient medium. Despite its 
importance, there is a factor of ~ 100 discrepancy between empirical and theoretical 
M of late-type O dwarfs, the so-called weak-wind problem. In this Letter, we propose 
a simple novel method to constrain M of runaway massive stars through observation 
of their bow shocks and Stromgrcn spheres, which might be of decisive importance 
for resolving the weak-wind problem. Using this method, we found that M of the 
well-known runaway 09.5 V star £ Oph is more than an order of magnitude higher 
than that derived from ultraviolet (UV) line-fitting (Marcolino et al. 2009) and is by a 
factor of 6 to 7 lower than those based on the theoretical recipe by Vink et al. (2000) 
and the Ha line (Mokiem et al. 2005). The discrepancy between M derived by our 
method and that based on UV lines would be even more severe if the stellar wind 
is clumpy. At the same time, our estimate of M agrees with that predicted by the 
moving reversing layer theory by Lucy (2010). 
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1 INTRODUCTION 

Mass-loss rate, M, is one of the key parameters defining the 
evolutionary sequence of massive stars and the end products 
of their evolution (e.g. Chiosi & Maeder 1986; Langer et al. 
1994). Although there is good agreement between the theo- 
retically predicted M of early-type O stars, i.e. stars more 
luminous than ss W 52 L@, and those derived empirically, 
there is a factor of 100 discrepancy between the theoretical 
and empirical M of late-type O dwarfs (e.g. Martins et al. 
2005b, 2012; Marcolino et al. 2009), the so-called weak-wind 
problem (see Hillier 2008 and Puis, Vink & Najarro 2008 for 
reviews) . 

M is also one of the main parameters defining the 
impact of massive stars on the interstellar medium (ISM) 
through the formation of stellar wind bubbles (e.g. Avedis- 
ova 1972; Castor, McCray & Weaver 1975) and bow shocks 
(Baranov, Krasnobaev & Kulikovskii 1971; Weaver et al. 
1977). The characteristic scale of these structures (typically 
tens of pc for bubbles and tenths of pc for bow shocks) is 
proportional to M and therefore, in principle, can be used 
to constrain this important parameter. The use of the wind 
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bubbles for this goal, however, is hampered because most 
massive stars reside in parent clusters so that their bub- 
bles are the result of the joint action of all massive stars in 
the cluster, which precludes us from constraining the wind 
parameters of individual stars. On the other hand, wind 
bubbles created by massive stars ejected into the field (the 
runaway stars) rapidly turn into bow shocks (Weaver et al. 
1977), whose sizes depend on individual characteristics of 
their underlying stars. This makes bow shocks an impor- 
tant tool for constraining M of massive stars (Gull & Sofia 
1979; Kobulnicky, Gilbert & Kiminki 2010) and thereby for 
resolving the weak-wind problem. 

The characteristic size of a bow shock depends on the 
stellar wind momentum rate (the product of M and the 
terminal velocity of the stellar wind, «oo), the space velocity 
of the star, and the number density of the ISM, n. The 
space velocities of nearby stars can be accurately determined 
by measuring their proper motions, parallaxes and radial 
velocities. Then Mvoo can be derived if n is known. 

In this Letter, we propose a simple novel method to 
constrain M of massive stars based on detection of well- 
defined H II regions around bow-shock-producing runaways. 
Detection of these structures around field O stars allows us 
to express Mvoo through observables (bow shock stand-off 
distance, Stromgren radius and space velocity) and the total 
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Figure 1. Left: SHASSA Ha image of the Hn region Sh2-27. The position of the ionizing star, £ Oph, is marked by a circle, while the 
direction of its peculiar (transverse) velocity is indicated by an arrow. The image is oriented with Galactic longitude (in units of degrees) 
increasing to the left and Galactic latitude increasing upward. Middle: IRAS 60 fim image of the same field with the bow shock generated 
by C Oph indicated by a dashed circle. The images were generated by the NASA's Sky View facility (McGlynn, Scollick & White 1998). 
Right: Spitzer 24 fim image of the bow shock around £ Oph. The orientation of the images is the same. At a distance of 112 pc, 1 degree 
corresponds to Rj 1.93 pc and 10 arcmin to £s 0.32 pc. 



Table 1. Summary of astrometric and kinematic data on f Oph (see text for details). 
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ionizing-photon luminosity, .S'(O), which is a well-established 
quantity for stars of known spectral type and luminosity 
class. We applied this method to the well-known runaway 
O star £ Oph, which according to Marcolino et al. (2009) 
belongs to the group of weak- wind stars (see Section[2] for a 
summary of relevant data on this star). The excellent con- 
sensus in the literature on Voo of £ Oph allowed us to esti- 
mate M of this star separately, which turns out to be inter- 
mediate between those based on Ha and ultraviolet (UV) 
lines (Section[3]). The discussion of this result and outlook 
are given in Section[4] 



2 ( OPH, ITS BOW SHOCK AND HII REGION 

C Oph (HD 149757, HIP 81377) is a nearby (« 112±| pc; van 
Leeuwen 2007), single (Garmany, Conti & Massey 1980), 
rapidly-rotating (vsini fa 400 kms -1 ; Howarth & Smith 
2001), runaway (Blaauw 1961), 09.5 Vnn (Morgan, Code 
& Whitford 1955; Lesh 1968) star with an impressive bow 
shock (Gull & Sofia 1979; van Buren & McCray 1988) im- 
mersed in the Hn region Sh2-27 (Sharpless 1959). 

In the left panel of Fig. [1] we present the Ha image of 
Sh2-27 originating from the Southern Hemispheric Ha Sky 
Survey Atlas (SHASSA; Gaustad et al. 2001), which shows 
an almost circular (ss 10° or 9.6 pc in diameter) Hn region. 
The middle panel of Fig.[T] shows the Infrared Astronomi- 
cal Satellite (IRAS) 60 (im image of the same field with an 
arcuate structure in the centre of the H II region, which is 
the bow shock generated by £ Oph (van Buren & McCray 



1988) . Using the archival Spitzer Space Telescope 24 fim im- 
age (Program Id.: 30088, PI: A. Noriega-Crespo) of the bow 
shock (see the right panel of Fig.[T} , we estimated the angu- 
lar separation between the apex of the bow shock and the 
star of « 5 arcmin, which corresponds to fa 0.16 pc. 

In Table[T] we give the parallactic distance, d, and 
the proper motion components, /x a cos5 and fig, of C, Oph 
(all from the new reduction of the Hipparcos data by van 
Leeuwen 2007), and the heliocentric radial velocity, « r ,hei) of 
this star (Wielen et al. 1999). To these data we added the 
components of the peculiar transverse velocity (in Galactic 
coordinates), Wj and Vb, the peculiar radial velocity, v r , and 
the total space velocity, v* = (vf + v% + v 2 ) 1 ^ 2 , of the star. 
To calculate v t , we used the Galactic constants Ro = 8.0 
kpc and 9 = 240 kms" 1 (Reid et al. 2009) and the so- 
lar peculiar motion {Uq,Vq,Wq) = (11.1, 12.2,7.3) kms -1 
(Schonrich, Binney & Dehnen 2010). For the error calcula- 
tion, only the errors of the proper motion measurement were 
considered. It follows from Tab le[T] that £ Oph is moving al- 
most in the plane of sky. 

The fundamental parameters of £ Oph were studied 
in numerous publications, of which the most recent are 
those by Mokiem et al. (2005) and Marcolino et al. (2009). 
There is excellent consensus on of the star: 1470 kms -1 
(Prinja, Barlow & Howarth 1990), 1550 kms -1 (Rcpolust, 
Puis & Herrero 2004) and 1500 kms -1 (Marcolino et al. 
2009). In what follows, we adopt — 1500 kms . Es- 
timates of M, however, are very different. The theoretical 
recipe by Vink, de Koter & Lamers (2000) predicts Myink ~ 
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1.29 x 10 7 M© yr 1 (for stellar parameters derived by Mar- 
colino et al. 2009), which is comparable to Mho; ~ 1-43 x 
10~ 7 M Q yr -1 , inferred by Mokiem et al. (2005) from synthe- 
sizing the Ha line using the FASTWIND code. A value two 
orders of magnitude lower, Muv ~ 1.58 x 10 -9 Mq yr -1 , was 
inferred by Marcolino et al. (2009) from CMFGEN model 
fits to the UV doublet of Civ AA1548, 1551, while Lucy 
(2010) predicted M LuC y ~ 1-30 x 10" 8 M yr -1 by using the 
updated moving reversing layer theory of Lucy & Solomon 
(1970). 

The fast rotation of £ Oph makes its spectral classifi- 
cation somewhat ambiguous. Although the widely accepted 
spectral type is 09.5 V, some other classifications were sug- 
gested as well. For instance, Conti & Leep (1974) classified 
this star as 09 V(e), while Herrero et al. (1992) prefer 09 III. 
Like Marcolino et al. (2009), we adopt the spectral type of 
09.5 V, so that 5(0) = 3.63 x 10 47 s~ 1 (Martins, Schaerer 
& Hillier 2005b). The basic parameters of £ Oph are sum- 
marized in Table[2] To this table we added our estimate of 
the mass- loss rate, M bs, based on the observed parameters 
of the bow shock and the H n region associated with £ Oph 
(see next Section). 



3 MASS-LOSS RATE OF < OPH 

The structure of an Hn region created by a moving star 
depends on the stellar velocity relative to the local ISM. 
For supersonically moving stars the number density within 
the H n region is comparable to that of the ambient ISM, 
while the radius of the Hn region is of the order of the 
Stromgren radius (e.g. Tenorio Tagle, Yorke & Bodenheimer 
1979), which is given by (e.g. Lequeux 2005) 



Rst 



( 35(0)^ 



1/3 



(1) 



where ob is the recombination coefficient of hydrogen to all 
but the ground state (S3 2.6 x 10 -13 cm 3 s _1 for fully ionized 
gas of temperature of 10 4 K). 

The stellar wind of a moving star interacts with the ISM 
and produces a bow shock ahead of the star. The minimum 
distance from the star at which the wind pressure is balanced 
by the ram and the thermal pressures of the ISM, the stand- 
off distance, is given by (e.g. Baranov et al. 1971) 
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where /i = 1.4 is the mean molecular weight, mu is the 
mass of a hydrogen atom, k is the Boltzmann constant, and 
T(~ 10 4 K) is the temperature of the (ionized) ISM. 

Eliminating n between equations JIJ and |2|, one has 
the 'observed' stellar wind momentum rate 
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where M — v t /c s is the isothermal Mach number and c s = 
(2kT / /imu) 1 ^ 2 is the isothermal sound speed (S3 10.9 kms -1 
for T = 10 4 K). For Rq = 0.16 pc, v, = 26.5 kms" 1 , 5(0) = 



3.63 x 10 47 s _1 , and Rst = 9.6 pc, one has from equation ([3| 
that M bs«oo ~ 2.1 x 10 26 gems" 2 . Thanks to the excellent 
consensus on of ( Oph (see Section[2}, one finds M of 
this star separately, M oba = 2.2 x 10~ 8 M Q yr _1 . This is 
S3 14 times larger than Muv, 7 times lower than Mhcx, and 
comparable to ML Ucy . The discrepancy between Muv and 
M b s would be even more severe if the wind is clumpy. For 
instance, Muv should be reduced by a factor of three (~ \fj) 
if the volume filling factor, /, in the wind is ~ 0.1 (Marcolino 
et al. 2009), so that M obs S3 40Muv- 

Using equations ((2)1 and (|3]), one finds n S3 3.6cm" 3 , 
which agrees well with estimates of the electron number 
density within Sh2-27 of ~ 3.0 cm" 3 (based on radio ob- 
servations, Ha surface brightness, and studies of interstellar 
absorption lines; Gull & Sofia 1979 and references therein) 
and the average line of sight electron and neutral hydro- 
gen densities towards £ Oph of S3 4.0 cm" 3 (Howk & Savage 
1999; scaled to d=112 pc). 

Note that M b s was derived for d — 112 pc, while Mar- 
colino et al. (2009) adopted a somewhat larger distance 
of 146 pc. An even larger value, d = 222 pc, was sug- 
gested by Megier et al. (2009), whose estimate is based on 
the empirical relationship between the strength of the in- 
terstellar Call lines and the distances to early-type stars. 
Since v* scales with the distance as d 6 (v* S3 31.2 and 
41.7 kms" 1 for d = 146 and 222 pc, respectively), one 
finds that M c b s tx d 17 , so that the larger the distance 
the larger the discrepancy between M b s and Muv, namely 
Mobs w 22 - 45Muv for d in the range from 146 to 222 pc. 
Correspondingly, Mn a remains larger than M b s by a factor 
of ss 2 to 4. Myink also depends on the adopted distance (via 
the stellar luminosity; Vink et al. 2000); this would decrease 
(increase) Myink by a factor of a few for d = 112 (222) pc. 

Similarly, it follows from equation (|3} that M b s oc 
5(0) 1/2 , so that M obs S3 21M UV (or S3 0.2Af HQ ) if < Oph is 
an 09 V star (Conti & Leep 1974), and M obs « 37Afuv (~ 
0.4M Ha ) if its spectral type is 09III (Herrero et al. 1992). 
Here we used the calibrations of Martins et al. (2005b). 

On the other hand, the possible leakage of ionizing pho- 
tons from the H II region (caused by porosity of the ISM; see 
Wood et al. 2005 and Sectionf?} would reduce our estimate 
of M bs- Conservatively assuming that a half of the photons 
escape from Sh2-27, one finds that M b s should be reduced 
by a factor of ~ 1.4. 

Note also that the possible existence of large-scale inter- 
nal flows in the H II region (e.g. caused by photoevaporation 
of molecular clumps in the ISM; see Section|3| also can af- 
fect the estimate of M G b s because Ro depends on the stellar 
velocity relative to the local ISM. Assuming that the char- 
acteristic velocity of the photoevaporation flows is of order 
of the sound speed in the ionized gas, i.e. ~ 10 kms" 1 , one 
finds that M bs might either be larger or smaller by a factor 
of up to S3 2 if the flow velocity is antiparallel or parallel to 
the vector of the space velocity of £ Oph. 

Finally, we note that the fast rotation of £ Oph might 
make the stellar wind anisotropic by increasing M in the 
equatorial zone (Friend & Abbott 1986) . This in turn might 
affect the geometry of the bow shock and the UV line and Ha 
diagnostic. The Ha diagnostic would also be unreliable if the 
Oe status of £ Oph is due to the presence of a circumstellar 
disk (see, however, Vink et al. 2009). The anisotropy effects 
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Table 2. Basic parameters of £ Oph (see text for details). 
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cannot be easily estimated, but are unlikely to change our 
result significantly. 

To summarize, the uncertainties in the distance and the 
spectral type might lead to the increase of M b s by a factor 
of ~ 1.5 to 3, while the leakage of ionizing photons from the 
Hn region and the presence of regular (e.g. photoevapora- 
tion) flows within the Hn region might reduce M b s by a 
factor of two. From this it follows that Af Q b s would remain 
at least an order of magnitude larger than Muv and might 
be comparable to or somewhat lower than Mua- 



4 DISCUSSION AND FURTHER WORK 

We estimated the mass- loss rate, M a bs, of the candidate 
weak-wind star £ Oph using the observed parameters of its 
bow shock and Hn region. We found that M b s is more 
than ten times larger than Muv inferred by Marcolino et 
al. (2009). (Recall that the difference between the two esti- 
mates would be even higher if the stellar wind is clumpy.) 
This finding supports the suggestion by Mokiem et al. (2007) 
that the use of only UV lines might significantly underesti- 
mate M of late-type O stars, because X-rays created by 
shocked wind regions can significantly change the wind ion- 
ization and thereby reduce the wind emission in the UV 
lines (see Martins et al. 2005a and Marcolino et al. 2009 
for more details). Moreover, the strengths of diagnostic UV 
lines could also be reduced (for a given M) if the wind is 
porous (clumpy) not only spatially, but also in velocity space 
(Sundqvist, Puis & Feldmeier 2010; Muijres et al. 2011). On 
the other hand, although the latter effect does not signif- 
icantly affect the Ha line, the Ha based estimates of M 
should be considered as upper limits because of the wind 
dumpiness (Repolust et al. 2004). 

The latter conclusion is in line with our finding that 
Mobs is « 7 times lower than Mho (derived by Mokiem et 
al. 2005 using the FASTWIND code). Note that the dis- 
crepancy between the two estimates could also be caused 
by the fact that the FASTWIND code tends to predict a 
stronger Ha absorption than the CMFGEN one (Puis et al. 
2005; Marcolino et al. 2009), which can affect the Ha based 
estimates of M. It is plausible therefore that the actual M 
of £ Oph is indeed somewhere between Muv and Mh q , as 
suggested by the work of Lucy (2010) and as found in this 
Letter. This inference along with the recent finding by Mui- 
jres et al. (2012) that the recipe by Vink et al. (2000) might 
over-predict M for late-type O dwarfs provides an avenue 
for resolving the weak-wind problem. 

Our estimate of M D b s was obtained under the assump- 
tion that £ Oph is running through a homogeneous ISM. The 
[Nil] /Ha and [Sii]/Ha line ratio maps of Sh2-27 (based on 
the Wisconsin Ha Mapper observations; Wood et al. 2005), 
suggest, however, that the ISM is porous, i.e. contains low- 
density voids. One might wonder therefore whether C, Oph is 



currently located within such a void so that its M is actually 
<C M ba and correspondingly comparable to Muv- 

To address this issue we note that the ionizing radiation 
of C, Oph could homogenize the ambient ISM by photoevap- 
oration of density inhomogeneities (e.g. Elmergreen 1976; 
McKee, van Buren & Lazareff 1984). Thus, our assumption 
of a homogeneous medium would be correct if the star com- 
pletely photoevaporates clumps within a region of radius 
^ P h ^ J?o on a time-scale, t p h, shorter than the crossing 
time of this region, i.e. t p h < R p h/v„. The photoevaporation 
time of a clump of mass M c \ located at a distance 7? p h is 
given by (McKee at al. 1984) 

X ( Me") (lie") (l0.9 kms- 1 ) ' (4) 

where n m is the mean density the local ISM would have if 
it were homogenized. The quite smooth appearance of the 
bow shock (whose transverse dimension is m 1 pc) implies 
that the local ISM is homogeneous on a scale of ~ 1 pc, so 
that our assumption would be correct if the mass of indi- 
vidual clumps is ^ 0.01 Mo (we assumed here that i? p h = 1 
pc and n m ~ n = 3.6cm- 3 ). It is difficult to make more 
quantitative arguments analytically, because this is a very 
non-linear process. We are currently using 3D radiation- 
magnetohydrodynamics simulations to investigate the ef- 
fects of a clumpy ISM in more detail (Mackey et al., in 
preparation). 

The current sample of Galactic candidate weak-wind 
stars (Martins et al. 2005, 2012; Marcolino et al. 2009) con- 
tains 22 stars (including £ Oph). We searched for bow shocks 
around all these stars using the Mid-Infrared All Sky Sur- 
vey carried out with the Wide-field Infrared Survey Explorer 
(WISE; Wright et al. 2010). This survey provides images 
in four wavebands centred at 3.4, 4.6, 12 and 22 /im, of 
which the 22 /an band is most suitable for detection of bow 
shocks (e.g. Gvaramadze et al. 2011; Peri et al. 2012). Be- 
sides £ Oph, the bow shocks were detected around four stars: 
HD 34078, HD 48099, HD 48279, and HD 216898. The first 
two stars are well-known runaways (see figs. 2 and 3 in Peri 
et al. 2012 for the WISE images of their bow shocks), while 
the runaway status of HD 48279 (suggested by its space ve- 
locity of ~ 30 kms- 1 and the presence of a bow shock) was 
recently established by Gvaramadze et al. (2012^0. The bow 
shock around HD 216898 is detected for the first time. 

Unfortunately, none of the four stars are associated with 



1 It is believed that HD 48279 is a member of the Mon OB2 as- 
sociation, while the orientation of the bow shock and the vector 
of the stellar peculiar transverse velocity imply that the star was 
injected into the association (cf . Gvaramadze & Bomans 2008) . 
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well-defined Hn regions. There are, however, several alter- 
native possibilities to estimate the local ISM density. For ex- 
ample, the pre-shock density could be derived through mod- 
elling (optical) spectra of the shocks (e.g. Dopita 1977) or 
from their emission measures, either in Ha or in the infrared 
(e.g., Bally et al. 2006; Kobulnicky et al. 2010). The latter 
possibility is especially attractive because bow shocks are 
detectable mostly in the infrared. Hydrodynamic modelling 
of bow shocks can also constrain the ISM density by com- 
paring observed instabilities, bow shock mass, and emissiv- 
ity to simulations (e.g. Mohamed, Mackey & Langer 2012). 
Further work in this direction is highly desirable in order to 
constrain M of a larger sample of 'weak-wind' stars. 

The paucity of bow-shock-producing 'weak-wind' stars 
with well-defined H n regions can be compensated by detec- 
tion of new examples of late O/early B type stars possess- 
ing these structures. Indeed, inspection of the SHASSA sur- 
vey revealed several circular Hn regions with central bow- 
shock-producing OB stars (Gvaramadze et al., in prepara- 
tion) , some of which (e.g. 8 Sco, r Sco) might belong to the 
group of the 'weak-wind' stars. Determination of M for these 
stars using different atmosphere models and their compari- 
son with Mobs (and M based on the work of Lucy 2010) is 
highly desirable as well. 

To conclude, detection of bow shocks (and well-defined 
H II regions) around field late- type O dwarfs provides unique 
possibilities for resolving the weak-wind problem, which in 
turn would have profound consequences for better under- 
standing the mass-loss mechanism of massive stars, their 
impact on the ambient ISM, and their evolution in general. 
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